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Abstract 



We report the observation of a narrow peak in the KsK~tt + invariant 
mass distribution in a sample of exclusive B — ► KKsK~ir + decays collected 
with the Belle detector at the KEKB asymmetric e + e _ collider. The measured 
mass of the peak is M = 3654±6(siat) ±8(syst) MeV/c 2 and we place a 90% 
confidence level upper limit on the width of F < 55 MeV/c 2 . The properties 
agree with heavy-quark potential model expectations for the n c (2S) meson, 
the n = 2 singlet S charmonium state. 
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Typeset using REVTpjX 



Major experimental issues for the charmed- quark anticharmed-quark (cc) charmonium 
particle system are the two cc states that are expected to be below open charm threshold but 
are still not well established: the radially excited n = 2 singlet S state, the r) c (2S) meson, and 
the n = 1 singlet P state, the h c (lP). The observation of these states and the determination 
of their masses would complete the below-threshold charmonium particle spectrum and 
provide useful information about the spin-spin part of the charmonium potential [|I|]. 

B meson decays provide an excellent opportunity for searching for the r] c (2S) and clari- 
fying its properties. They are a copious ^(I'S') source: the decays B — > Krj c {lS) have been 
observed by CLEO |2J, BaBar [Q, and Belle 0] with relatively large branching fractions: 
B{B — > Kr] c (lS)) ~ B(B — > KJ/ip) ~ 1 x 10~ 3 . (In the following, we use r] c to denote 
the 77 C (1S').) Moreover, in the case of the triplet charmonium S states, B meson decays to 
the radially excited ip{2S) are nearly as common as those to the n = 1 J/ip radial ground 
state: B(B + K+ip(2S))/B(B + -> K+J/^p) ~ 0.6 @. Thus, it is reasonable to expect the 
decays B —>■ Krj c (2S) to occur at a rate comparable to that for B — > Kr] c . Unlike the J/i/j 
and ip(2S), where hadronic decays proceed via highly suppressed three-gluon intermediate 
states, the r\ c and r) c (2S) decay via less-suppressed two-gluon processes. As a result, inter- 
charmonium transitions are not very important and the hadronic decay branching fractions 
for the two states are expected to be similar ||. Thus, any final state that shows a strong 
B — > Kr] c signal is a promising channel for an r] c (2S) search. 

A simple application of heavy-quark potential models [0] predicts a ip{2S)-rj c {2S) mass 
splitting that is smaller than that for the ground-state J / ip-r] c splitting because of the smaller 
value of the wave function at zero cc separations and the running of the QCD coupling 
strength a s (M 2 ). These models predict an T] C (2S) mass in the range 3625 < M Vc ^s) < 
3645 MeV/c 2 . Similar factors result in the expectation that the r] c (2S) total width is some- 
what smaller than that of the r\ c . 

The Crystal Ball group || reported an excess of E 1 ~ 91 MeV gamma rays from inclusive 
ip(2S) — > 7 A decays, and interpreted this as possible evidence for the t] c {2S) with mass 
3594 ± 5 MeV/c 2 . This result implies a tp(2S)-r] c (2S) mass splitting that is considerably 
larger than heavy-quark potential model expectations. The result has not been confirmed 
by other experiments 0. 

In this letter we report a search for the r] c (2S) produced via the processes B + — > 
K + r] c (2S) and B° -> K s r] c (2S), where rj c (2S) -> K s R-n + 0. We concentrate on this 
final state because it is a strong decay channel for the r] c (B ~ 1.8%), has low combinatorial 
backgrounds, and, since the final state contains all charged particles, is reconstructed with 
good resolution. Moreover, the process ip{2S) — > KsK~n + is strongly suppressed, and the 
background in this channel from B — > Kip(2S) decays is expected to be less than 0.1 events. 

The search uses a 41.8 fb _1 data sample collected with the Belle detector [0] at the 



KEKB e+e~ collider pf operating at the T(45) resonance (y/s = 10.58 GeV). The data 



sample contains 44.8 million BB meson pairs. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector, a 50-layer central drift chamber (CDC), an array of aerogel 
threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation 
counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals located 
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inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect Ki mesons and to identify muons. 



The detector is described in detail elsewhere 111 



We select events with K + KsK ± tt t or KsKsK ± n T combinations. Here the charged 
kaon (pion) tracks are required to originate from within 5r < 0.3 cm and \5z\ < 2.2 cm of 
the run-by-run determined interaction point (IP) in the transverse (r<p) and beamline (z) 
directions, respectively. In addition, they must be positively identified as kaons (pions) by 
the combined information from the ACC, TOF and CDC dE/dx measurement. Candidate 
Ks — > 7r + 7r~ decays correspond to pairs of oppositely charged tracks with invariant mass 
within 12 MeV/c 2 (3<r) of M K o that originate from a common vertex that is displaced by 
more than 0.3 cm from the IP. The direction of the Ks momentum vector is required to be 
within 0.2 rad of the direction between the IP and the position of the displaced vertex. 

Candidate B mesons are reconstructed using the energy difference AE = E B ms — E^ m 



and the beam-energy constrained mass M^ c = J (E^ m ) 2 — (pg 118 ) 2 , where E^ m is the 
center of mass (cms) beam energy, and E B ms and p c ™ are the cms energy and momentum of 
the B candidate. The signal region is defined as 5.271 < M bc < 5.287 GeV/c 2 and \AE\ < 
0.040 GeV, which correspond to ±3cr from the central values for both quantities. 

In order to suppress background from the e + e~ — > qq continuum (q = u, d, s & c), 
we form a likelihood ratio from two variables. One is a Fisher discriminant determined 
from five modified Fox- Wolfram moments |13[], the cosine of the angle formed by the thrust 
axis of the candidate B — > KKsK~tt + tracks and that of the remaining tracks in the 
event, and the sum of the absolute values of transverse momenta of particles relative to the 
B candidate's thrust axis with angle larger than 60° normalized by the sum of the total 
momenta. The coefficients of the Fisher discriminant are chosen to optimize the separation 



between signal and continuum Monte Carlo (MC) events [14]. The other is the cosine of 
the angle between the B candidate flight direction and the beam axis in the T(45) rest 
frame (cos#b). Normalized probability density functions (pdfs) formed from the Fisher 
discriminants and the cos 9b distribution are multiplied to form likelihood functions for 
the signal {C sig ) and continuum (C cont ) processes. We select events with a likelihood ratio 
LR = C si g/ (C si g+C cont ) > 0.6, which was determined by optimizing S/ V S + B (S and B are 
signal and background, respectively) for Monte Carlo simulations of the process B — > Kr] c , 
where t] c KsK~tt + . 

We reduce potential backgrounds from B — > D(D S )X decays by rejecting D and D s 
mesons with the requirements \M K7r -M D \ > 10 MeV/c 2 and \M KsK +-M Ds \ > 10 MeV/c 2 . 
The decay rj c (nS) — > K*(890)K is suppressed by an angular momentum barrier; in order to 
reduce backgrounds from other B meson decay modes with minimal loss in signal, we reject 
events with a K* candidate with the requirement \M Kn — M K *\ > 50 MeV/c 2 . 

Figure [I] shows the M^ c projections of events in the \ AE\ < 0.040 GeV signal region for 
twenty five Mx s Kn mass bins, each 40 MeV/c 2 -wide and with central values ranging from 
2840 through 3800 MeV/c 2 . The mass bins of Figs. |I|(d) and (e) straddle M Vc and clear 
peaks corresponding to B — » Kr) c , r\ c — * K s K~ti + decays are apparent. Figures [l|(u) and 
(v), which cover a region near the expected mass of the t] c (2S), also show distinct B meson 
signals. 
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FIG. 1. The Mb c projections for 40 MeV/c 2 bins of Mk s k-ki with central values ranging from 
2840 (a) to 3800 MeV/c 2 (y). Only events with |AE| < 40 MeV are included; the charged and 
neutral B decay modes are combined. 

We perform simultaneous fits to each of the Mb c distributions of Fig. || and the corre- 
sponding AE distributions for events in the M^ c signal region (not shown). The fits use 
Gaussian functions with MC-determined widths to represent the signals; the areas of the 
Mbc and AE signal functions are constrained to be equal. The Mb c background is modeled 
by a smooth function that behaves like phase-space near the kinematic end point 15|| ; for 
the AE background, we use a second-order polynomial. As an example, the results of the fit 
to the Mbc and AE distributions of the Mx s Kn — 3640 MeV/c 2 bin are shown in Figs. 0(a) 
and (b), respectively. 



6 




3620<M(K s Ktt)<3660 MeV/c 2 

15 



10 



5.300 





-0.20 



- 1>) 

— "I 


- 

In 'iLi 







0.20 



M 



0.00 

,. bc (GeV/c*) AE (GeV) 

The (a) M hc and (b) AE projections for the M Ks Km = 3640 MeV/c 2 mass bin. The 
curves are the results of the simultaneous fit described in the text. 



FIG. 2. 



The signal yields extracted from the simultaneous fits to the different KsK~tt + mass 
bins are plotted vs. Mk s k-k in Fig- H, where, in addition to a prominent r\ c peak and a hint 
of a J/^, a clear peak at higher mass is evident. We identify this as a candidate for the 
f] c (2S). Between the peaks is a non-zero, non-resonant contribution. The curve in Fig. ^ is 
the result of a fit with simple Breit-Wigner functions that represent the rj c and candidate 
r] c (2S), a Gaussian function with mass and width fixed at the J/ip values, and a second-order 
polynomial to represent the non-resonant contribution. These functions are convolved with 
a Gaussian resolution function with a MC-determined width of a = 15 MeV/c 2 . 
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FIG. 3. The distribution of signal events from the simultaneous fits to M\> c and AE for each 
KsKtt mass bin. The curve is the result of the fit described in the text. 

The fit values for the event yields, masses and total widths of the rj c and the T] C (2S) 
candidate state are listed with their statistical errors in Table |. The fit value for the r\ c 
mass is in good agreement with the world-average value of M Vc = 2979.8 ± 1.8 MeV/c 2 ||; 
the value for the r) c width is consistent, within its rather large errors, both with the existing 

and the recent CLEO result of 26 ±6 MeV/c 2 



13.2±1;1 MeV/c 2 
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world average of r* ot 

° '/c 

The sum of the observed events in the three mass bins in the signal region (i.e., centered 
around M(KsKtt) = 3640 MeV/c 2 ) is 56, while the integral of the second-order polynomial 
over the same interval gives a non-resonant expectation of 21 ± 2 events. The probability 
for this to fluctuate up to 56 events is ~ 10~ 8 , which corresponds to a signal significance of 
more than 6a. 

The fitted mass of the candidate r) c (2S) is substantially above the Crystal Ball mass 
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value and consistent, within errors, with the upper end of potential model expectations. The 
systematic error on the mass is evaluated by redoing the analysis using different likelihood 
ratio selection requirements, 50 MeV/c 2 -wide bins, bins with central values shifted by half a 
bin-width, and with different values of the experimental resolution. The maximum change 
in the fitted mass value is 8 MeV/c 2 , which is taken as the systematic error. The limited 
statistics and the resolution precludes a precise measurement of the width. However, we can 
establish a 90% confidence level upper limit of F < 55 MeV/c 2 . 

Monte Carlo simulations indicate that the acceptance is very nearly constant for the 
KsK~n + mass region covered by this measurement JT7J. Thus, the ratio of product branch- 
ing fractions for the r\ c and rj c (2S) is just the ratio of event yields: 

B(B : IB K ' U( 1 S ^Vrf^ = ^ ± 0-12 ± 0.05. (1) 
B(B -> Jf!fc)B(!fc - K s K-k+) y ' 

where the first error is statistical and the second systematic. The systematic error is de- 
termined from changes in the ratio observed for different binning, values of resolution, and 
functions used to model the non-resonant contribution. 

In summary, we observe a peak in the KsK~n + mass from exclusive B + — > K + KsK~tt + 
and B° — > KsKsK~tt + decays with mass and width values 

M = 3654 ± 6 ± 8 MeV/c 2 
T < 55 MeV/c 2 ; 

these are consistent with expectations for the B — > Kr] c (2S), where r] c (2S) — > KsK~tt + . In 
addition, the product branching fraction is comparable in magnitude to that for the tj c , also 
in agreement with expectations for the rj c (2S). The observed properties of this system lead 
us to conclude that we have observed the r) c (2S). 
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Science, and Technology of Japan and the Japan Society for the Promotion of Science; 
the Australian Research Council and the Australian Department of Industry, Science and 
Resources; the National Science Foundation of China under contract No. 10175071; the 
Department of Science and Technology of India; the BK21 program of the Ministry of 
Education of Korea and the CHEP SRC program of the Korea Science and Engineering 
Foundation; the Polish State Committee for Scientific Research under contract No. 2P03B 
17017; the Ministry of Science and Technology of the Russian Federation; the Ministry of 
Education, Science and Sport of the Republic of Slovenia; the National Science Council and 
the Ministry of Education of Taiwan; and the U.S. Department of Energy. 
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TABLES 



TABLE I. Results of the fit to the data points in Fig. |3|. Only statistical errors are listed. 



Peak 


Events 


Mass (MeV/c 2 ) 


L tot (MeV/c 2 ) 




104 ± 14 


2979 ± 2 


11 ± 11 


Vc(2S) 


39 ±11 


3654 ± 6 


15+ 24 

i0 -15 
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